Stability of the Smectic Quantum Hall State: A Quantitative Study 
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We present an effective elastic theory which quantitatively describes the stripe phase of the two- 
dimensional electron gas in high Landau levels (N > 2). The dynamical matrix is obtained with re- 
markably high precision from the density-density correlation function in the time-dependent Hartree- 
Fock approximation. A renormalization group analysis shows that at T = 0, as the partial filling 
factor Av = v — \y\ moves away from 1/2, the anisotropic conducting state may undergo quan- 
tum phase transitions: stripes may get pinned along their conducting direction by disorder, or may 
lock into one another to form a two-dimensional crystal. The model predicts values of Av for each 
transition. The transitions should be reflected in the temperature dependence of the dissipative 
conductivity. 

PACS numbers: 73.40.Hm, 73.20.Dx, 73.20.Mf, 73.40.Kp 



In a strong magnetic field perpendicular to the plane of 
a two-dimensional (2D) electron or hole system, the en- 
ergy spectrum is characterized by discrete Landau levels 
(LLs) separated by the cyclotron energy %u> c = TieB/mc 
and the Zeeman energy g^isB. Since the degeneracy of 
each LL, A/^, is proportional to the magnetic field, for 
sufficiently strong magnetic fields, only a small number 
of low-lying LLs are occupied. In this situation, the ki- 
netic energy is practically quenched, sometimes leading 
to interesting strongly correlated liquid ground states |Q . 
For example, it is now very well established that the 
ground state of a sufficiently clean system is given by 
the fractional quantum Hall fluid |^| if the filling factor 
v = N (AT = number of electrons or holes), which 
measures how many LLs are filled, is close to some ratio- 
nal numbers such as 1/3, 1/5, 2/3, 2/5, etc. 

However, recent experiments on high- mobility 2D elec- 
tron H and hole [Q systems have revealed that the trans- 
port properties are qualitatively different for higher fill- 
ing factors (y > 4 for electrons and v > 2 for holes.) 
One of the most remarkable findings is that the longi- 
tudinal resistivities are highly anisotropic near half inte- 
ger filling factors. As a natural explanation, it has been 
suggested that the electrons form unidirectional charge 
density waves or "stripes" , predicted earlier by Koulakov 
et al. H and Moessner and Chalker An important 
theoretical development M was the observation that the 
states of this system may be classified by their symme- 
tries, which are highly analogous to those of liquid crys- 
tals. The possible states include stripe crystals, smec- 
tic, and nematic phases |^|. Of these, the smectic shows 
promise of explaining the anisotropic transport data at 
very low temperature |9| . 

While some progress has been made in understanding 
the temperature dependence of transport in this system 
pfj| , the true groundstate when quantum fluctuations 
are included remains a subject of debate. It has been 
shown |ll|,^,^2| that continuous quantum phase transi- 
tions may occur among different possible states when the 
system parameters are varied. To assess whether this 



happens requires a knowledge of the parameters enter- 
ing the effective theory from microscopic calculations. In 
one such study, MacDonald and Fisher, using an elas- 
tic edge state model, found that both interstripe locking 
interactions p| and pinning by weak random disorder 
are relevant perturbations at all Av |l3|. However, they 
argued that both perturbations are extremely small in 
any experimentally accessible systems, and computed the 
anisotropic resistivities using a semiclassical Boltzmann 
transport theory. 
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FIG. 1. The anisotropic lattice used in the effective elastic 
model. 

In this work, we focus on filling factors away from 
Av =1/2. In recent work @, it was found that in the 
Hartree-Fock approximation, uniform stripe states are 
unstable to the formation of density modulations along 
the stripes — i.e., to the formation of a stripe crystal. 
Interestingly, the collective modes 0] of the modulated 
stripe state obtained from a time-dependent Hartree- 
Fock approximation (TDHFA) |L5| shows that the mo- 
tion of density deformations in the low energy modes 
]l6| looks strikingly similar to that of phonon modes in a 
highly anisotropic lattice (see Fig. Q). Motivated by this 
observation, we develop an effective elastic lattice model 
which provides a quantitative description of the low tem- 
perature behavior of the stripe phase. As described be- 
low, from the density-density correlation function in the 
TDHFA, we may numerically obtain a dynamical matrix 
which reproduces the low-energy microscopic behavior of 
the system with remarkably high precision. Using the 
result, we perform a renormalization group (RG) analy- 
sis on two important perturbations: interactions among 
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stripe modulations, which presumably stabilize a stripe 
crystal, and disorder, which may pin the electrons and 
render the smectic insulating. In contrast to uniform 
stripes we find for modulated stripes that, both these 
perturbations are irrelevant close to Av — 1/2, and that 
they become relevant at different critical values as Av 
moves away from 1 /2. Either of these Kostcrlitz-Thouless 
(KT) transitions represent a metal-insulator transition, 
and could be identified experimentally by measuring the 
filling factor dependence of the activation gap on the in- 
sulating (i.e., quantized Hall effect) side of the transition. 

In our simple lattice model, we assume that each lattice 
site R is occupied by an object with a local density profile 
or a "form factor" /(r), and that the only dynamic vari- 
able is the displacement u(R) of the object center. The 
clastic energy associated with the displacements may be 
described by a Hamiltonian in Fourier space as 

Ho = ^^2^2ui(q)D a/3 {q)uf 3 {q), (a,(3 = x,y). (1) 

q a/3 

Since the Hamiltonian is Hermitian and invariant under 
a rotation by 7r, all dynamical density matrix elements 
are real and D xy (q) = D yx (q), leaving us only three 
independent real parameters per wave vector. Below, we 
will treat D a p as a fitting parameter, chosen to match the 
density-density correlation function |l7]] obtained from 
the TDHFA of the underlying microscopic fermion model 
in the partially filled uppermost LL. 

Since the density-density correlation function is com- 
puted at the one-phonon level in the TDHFA 
a proper matching procedure requires the corresponding 
quantity in the lattice theory to be defined also within 
a harmonic approximation. More specifically, assuming 
q • u is small, the density fluctuation operator in the lat- 
tice model may be written as |18 



Sn(q) ee n(q) - (ri(q)) 



/(q) 



q • u(q) +Yj [ dr e tG -[ r - u Ml- lq - r 



, (2) 



where G is a reciprocal lattice vector and (• ■ •) denotes an 
average over the above elastic Hamiltonian. Ignoring the 
higher-order second term, the density-density correlation 
function may be written as 

X (q, q', t) ee (TSn(q, r)Sn(-q',0)) 
/(q)/*(q') 



{Nabf 



^ Q ^(T MQ (q,T)4(q',0)), (3) 

a/3 



where T is the time ordering operator and a and b are 
the lattice constants as in Fig. |l]. Due to the lattice 
symmetry, u a (q) is defined only within the first Brillouin 
zone, and the above expression vanishes unless q' = q+G 
for any reciprocal lattice vector G. 

Using the fact that the single-LL-projected displace- 
ment operators obey [u x (q), u y (q')] = — i^<5 q + q <,o {Ib — 



^Jhc/eB is the magnetic length), the above correlation 
function is easily computed, and may be written in terms 
of the Matsubara frequency ioj n as 



X(q+G,q+G',iw„) = 

VF(q;G,G') IF*(q;G,G') 



(4) 



lUJ n — LUq iuj n + Wq 

where W(q; G, G') are weights that depend on D a p, and 



i ~ ^b\ D xx D yy — D% is the eigenmode frequency. 
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FIG. 2. The phonon eigen frequency for Av = 0.4 at sev- 
eral wave vectors. The solid line is the result from the TDHFA 
and the circles are the result of fitting using the lattice model. 

For each wave vector q in the first Brillouin zone, we 
need to compute x f° r a t least three reciprocal lattice 
vectors to fit all unknown parameters in the problem. 
Assuming the stripes are along the s-axis as in Fig. [I], 
we use |9| G,G' = 0,±G (G = y27r/a), to find the 
three dynamical matrix elements D xx , D xy , and D yy , as 
well as the three form factors /(q), and /(q±Go). (Due 
to the inversion symmetry, the /'s are real.) For con- 
creteness, we set \v\ = 6 (the third lowest LL) in our 
calculations. In our matching procedure, we focus on the 
3x3 Hermitian matrix W(q; G, G') with G, G' = 0, ±G 
provided to us by the TDHFA. Together with the mode 
frequency w q , this gives ten parameters, of which we use 
six to produce the dynamical matrix and form factors. 
(In practice, it is easiest to use the off-diagonal elements 
of the Hermitian matrix.) The remaining four are used 
to estimate an errorbar and to examine the validity of 
the fitting. For example, Fig. || shows cj q obtained from 
the TDHFA and the lattice theory. For all the quanti- 
ties, the agreement is quite impressive. The relative error 
is typically ~ 10 -8 unless w q is very small, and always 
remains below ~ 10 . 

The numerical solution of D a p agrees remarkably well 
with a harmonic theory of a 2D charged smectic system 
with long-range Coulomb interaction such as discussed 
in Ref. [Q . Specifically, the functional forms of the dy- 
namical matrix elements may be written for small q as 
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D x 
D, 
D 



llll 



ll + ill (I, 
QxQy + QxQy/q, 

Qy + ql/q + qi, 



(5) 



where the terms proportional to l/q are from the long- 
range Coulomb interaction, the q x term from bending 
energy of each stripe, and all the rest from effective short- 
range interactions. 

Note that w q = at q x = 0. In fact, the only true gap- 
less mode of the stripe crystal is at q = 0. However, the 
gaps are extremely small if q x = 0, for any q y 14 . (In 



these modes, the stripes slide with respect to one another 
fL6[|.) We are thus able to model the results of the TD- 
HFA as a charged smectic, despite the formal presence of 
gaps away from q = 0. Obviously, these nearly gapless 
modes greatly affect the low temperature behavior; the 
RG exponents depend on all g y 's up to the first Brillouin 
zone boundary. Thus, a reliable elastic Hamiltonian must 
match the microscopic system at short wavelengths per- 
pendicular to the stripes. Along the stripes, however, 
only the long- wavelength behavior is important. 




FIG. 3. The renormalization group exponents for locking 
(2 — x K /2) and pinning along the stripes (3/2 — iv). The ex- 
ponents monotonically decreases as Av approaches 1/2. Er- 
rorbars are too small to draw in scale. 

Using our quantitative elastic theory, we examine some 
important perturbations and their effect at low tempera- 
tures by computing the RG exponents. We first consider 
locking between stripes, which controls the transition be- 
tween the smectic and the stripe crystal phases. The 
perturbative action may be written as |0] 



Sx = 



-A J dxdr cos — [u x (x, j, t) 



i x {x,j - 1,t)] . (6) 



Here we represent the displacements in real space, with 
the j labeling the stripes. For convenience, we have 
adopted a continuum approximation in the s-axis, which 



should not affect the small q x behavior. A standard RG 
analysis Efl] yields the flow equation 



where 



2al% 



dq y (l 



cos q v a) lim 



Dyy(q)q x 



(7) 



(8) 



The exponent is plotted as a function of Av in Fig. [|. 
Note that if Av > 0.38, then 2 — x K /2 becomes negative 
and locking is irrelevant 

Another important perturbation is that of weak dis- 
order, which may pin the stripes. We model this with, 
perturbative action for white noise disorder, 



S v = / drdTV(r)n(r), 



(9) 



where the disorder average is given by V(r)V(r') = 
V 2 abS(r - r'). Using the density in Eq. (g), the RG 
analysis is straightforward. If the G in the second term 
of Eq. (g) is parallel to the stripes, the flow equation for 
the most relevant operator G x — 2ir/b is given by 



di 



V . 



with 



< f 'a 



lim — — 

o UJ, 



(10) 



(ii) 



As plotted in Fig. S, if Av > 0.43, the exponent is nega- 
tive and pinning along the stripes is irrelevant. 

The RG analysis is more complicated if G is perpen- 
dicular to the stripes. We find, for any G = (0, G y ), that 
the free energy F = — In / X>uexp(— So — 5V) computed 
using Eq. (||) diverges at low temperatures as J 1 " 2 / 5 for 
any Av. This indicates that pinning across the stripes 
is always relevant. Our interpretation of this is that the 
stripes will be trapped in channels; however, they are still 
free to move along the channels so that this would not 
spoil the phase transitions we found above. 

Using particle-hole symmetry for Av > 1/2, the above 
results imply that for 0.43 < Av < 0.57, the stripes 
are pinned only in the the direction perpendicular to 
the stripes at T = 0, supporting the existence of an 
anisotropic metallic state This is the main re- 

sult of this Letter, and is in qualitative agreement with 
experiment [[|. 

Our results suggest that the anisotropic transport 
properties observed in experiment may well represent a 
new and unusual metallic state of the quantum Hall sys- 
tem, separated from insulating (quantized Hall) states by 
quantum phase transitions as a function of v, on either 
side of Av = 1/2. One direct probe of this is the ac- 
tivation energy A of the diagonal transport coefficients 
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Pxx, Pyy in the quantized Hall state. As the transition is 
approached from below, the gap is controlled by the grow- 
ing correlation length £, which for a KT transition has 
the characteristic form A ~ ~ exp {— C/^J\v c — v\}, 
where C is a non- universal constant. The observation 
of such behavior in finite temperature studies away from 
Av = 1/2 would constitute direct evidence of a new type 
of phase transition for these systems. 

We conclude with a few final remarks. Our simple 
model was analyzed only in the third lowest LL and it 
included neither finite thickness of the 2D layer nor more 
complex effects such as spin-orbit coupling. We expect 
that the qualitative picture presented will be robust, al- 
though the precise location of the phase transitions will 
surely be affected. Although in our specific model, pin- 
ning preempted locking, the order of transitions may well 
be reversed when these effects are included. We note also 
that the method we have developed here works very well 
away from Av = 1/2, but breaks down close to this value. 
The breakdown occurs for \Av — 1/2 < 0.01, and is re- 
flected in rapidly growing error bars for our matching 
procedure in that small range. This is in part caused by 
the charge motion of the collective modes found in TD- 
HFA becoming more complicated near half-filling, tak- 
ing on a mixed character of both edge states and lattice 
phonons. Apparently our simple elastic model does not 
capture this complicated behavior. 

In summary, we have developed a method for generat- 
ing a quantitative elastic model of quantum Hall stripes 
from microscopic TDHFA calculations. An RG analysis 
shows that an anisotropic conducting (smectic) state may 
be stable against crystallization and pinning near filling 
factor 1/2, and may be destabilized by either of these 
away from this filling in a continuous quantum phase 
transition. 
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